The convex hull of a regular set of integer vec- 
tors is polyhedral and effectively computable 
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Abstract 

Number Decision Diagrams (NDD) provide a natural finite symbolic representation 
for regular set of integer vectors encoded as strings of digit vectors (least or most 
significant digit first). The convex hull of the set of vectors represented by a NDD 
is proved to be an effectively computable convex polyhedron. 
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Presburger arithmetic [Pre29l ] is a decida ble logic used in a large range of ap- 
plications. Different techniques GBD02I ] and tools have been developed for 
manipulating the Presburger- definable sets (the sets of integer vectors satis- 
fying a Presburger formula ): by working directly on the Presburger-formulas 
(implemented in Omeg a |Ome| ) . by using semi- linear sets GS66] (imple- 
mented in Brain RV02| ) . or by using NDD (automata that represent regular 
sets of integer vecto rs enc o ded as strings of digit vectors, least or most sig - 
nifican t dig it first) |Boi98l . |WB95l . |BC96| (implemented in Fast |BFLP03 



Lash Las| and CSL-ALV|BB03]). Presburger-formulas and semi-linear sets 
lack canonicity: there does not exist a natural way to canonically represent 
a set. As a direct consequence, a set that possesses a simple representation 
could unfortunately be represented in an unduly complicated way. Moreover, 
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deciding if a given vector of integers is in a given set, is at least NP-hard 



Ber77l . |GS66| . On the other hand, a minimization procedure for automata 
provides a canonical representation for NDD- definable sets (a set represented 
by a NDD). That means, the NDD that represents a given set only depends 
on the set and not on the way we have computed it. For this reason, NDD are 
well adapted for applications that require a lot of Boolean manipulations like 
model-checking. 

Verification of systems with unbounded integer variables is undecidable in 
general. That explains why we are interested in over- approximating the reach- 
ability set of such a system. By computing the convex hull of the set of in itial 
states of such a system and by using a widening operator |CH78 . HPR97 



an 

over-approximation of the set of reachability set can be effectively computed. 

In this presentation, the convex hull of a set of integer vectors represented 
by a NDD is proved to be a convex polyhedron. That shows that it can be 
finitely represented as a finite intersection of half-spaces or dually as a finite 
set of rays. Indeed, we provide an exponential time algorithm that effectively 
computes this convex hull (the exact complexity remains open). 

This result is obtained by first proving that "the convex hull" of the language 
cr*.<T2 is equal to the convex hull of o~\.o\ for any pair of words (<7i, a%). From 
this commutativity result, we deduce that the convex hull of any regular lan- 
guage L, is equal to the convex hull of a finite union of regular languages of 
the form u>o.cr* . . . w n -i.a^.w n . 



1 Closed sets and convex sets 



Recall that the scalar product of two real vectors x, y G M. m where m > 1 is the 
real (x, y) = Y^i x [i]-v[i] where x[i] G M corresponds to the ith component 
of x. We denote by \x\ 2 the norm |x| 2 = J {x, x). The open ball centered in 
x G M m with a radius e > is the subset B Xfi = {y G W m ; \x — y\ 2 < e}. Recall 
that a subset X C M m is said open if for any x G X there exists e > such 
that B Xie CI. A closed set X is a subset of M m such that difference R m \X is 
open. Recall that any infinite or finite intersection of closed sets is closed and 
any subset X is included into a minimal (for the inclusion) closed set, called 
the closure of X. We denote by cl : P(R m ) — > V(M. m ) the function such that 
cl(X) is the closure of X for any X C R m . 

An half-space H is a subset of real vectors W 71 such that there exists a G M m 
and ceM satisfying H = {x G W 71 ; (a, x) + c#0} where # G {>, >}. Recall 
that such an half space H is closed if # is equal to > and it is open if # is 
equal to >. 
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We denote by R+ and K_ respectively the set of non-negative reals IR + = {iG 
K; x > 0} and the set of non-positive reals R_{x G M; x < 0}. 

A convex set is a finite or infinite intersection of half-spaces. The convex hull 
of a subset X C ]R m is the least (for the inclusion C) convex set that contains 
X. We denote by conv : V(R. m ) — > P(R m ) the function such that conv(X) is 
the convex hull of X for any X C M m . Recall that a vector y is in conv(X) 
if and only if there exists a finite sequence (xi)x<i< n of n > 1 vectors in X 
and a sequence (*i)i<*<n of n reals in IR + such that Z)r=i = 1 an d such that 
y = Yh=i U- X i- Recall that the closure of a convex set remains a convex set. 

A convex set C is said polyhedral if C is equal to a finite intersection of closed 
half-spaces (in particular, a polyhedral convex set is closed). Recall that any 
polyhedral convex set P can be represented by a finite set of rays R C R m x R + 
such that P = P(P) = {x G M m ; (x, 1) G C(P)} where C(P) C R w x R + is 
the polyhedral cone defined by the following equality: 

C{R) = {Y^U-r; t r G M+} 

re-R 

Recall that for any pair (Pi, P2) of polyhedral convex sets respectively repre- 
sented by a pair of finite set of rays (Pi, R%), the convex set cl o conv(Pi U P2) 
is polyhedral and represented by the set of rays Pi U P2. 



2 Regular sets of integer vectors 



Let us consider an integer r > 2 called the basis of the decomposition and an 
integer m > 1 called the dimension of the represented vectors. A digit vector b 
is an element of the finite alphabet S r m = {0 . . . r — l} m . The vector p(a) G N m 
associated with a word a = b\ . . . b n of n > 1 digit vectors 6j G S r ™ is defined 
by p{o~) = Yh=i f l ~ X -bi- We naturally define p(e) = (0, . . . , 0), also written 0. 

The set X represented by a language L C E* m is defined by X = p(L) = 
{p(cr); a G L}. If L is regular (that means accepted by a /imte automaton) , 
the set X is naturally said regular. Let us recall that r egular s ets of v ectors can 
be efficiently manipulated with finite automata (see and they 



correspond to the sets defined by a formula in the first order logic (N, +, <, K) 
where K is the valuation function in base r defined by y = V r (x) if and only 



if y is the greatest power of r that divides x [BHMV94I . 



Example 1 Consider the following automaton A + with basis r = 2 and di- 
mension m = 3 depicted below. Intuitively, this automaton represents the set 
of vectors (x, y, z) G N 3 such that x + y = z where the state qi corresponds to 
the carry i G {0, 1} of the addition. 
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3 The convex hull of a regular set of integer vectors 

The main result of this paper is proved in this section. We show that the 
closure of the convex hull of a regular set of integer vectors is polyhedral and 
represented by a set of rays effectively computable in exponential time from 
any regular expression that defines this regular set. 

As p(aw) = r' CT '.p(to) + p(cr) for any pair of words (a, w), we introduce the 
function T a : R m — > M. m defined by T a (x) = r'°"'.x + p(cr). Remark that any 
regular language L can be decomposed into a finite union of regular languages 
of the form o n +\.L* n .o n . . . L\.o~\ where cr; G S* m and Li C £* m . The following 
lemma is a first step toward the computation of cl o conv o p[L). 

Proposition 2 Let us consider a language L = cr„ +1 .L* .cr n . . . L\.o\ where 
n > 0, Oi G S* m and Lj C S* m . We have the following equality: 

cloconvop(L) = r CTn+ ,.. CT1 oclo^{(0, . . . ,0)} U \J R-T" 1 CTi o cl o conv o 
where a — > ^(a) is partially defined over T^ m , by the following equality: 



PROOF. We denote by C(a n+1 , L n , cr n , . . . , L 1; ai), the following set: 

<T„-|_1 ...<T1 

o cl o {(0, . . . , 0)} U (J R-T-l ai o cl o conv o £(L*) 
V i=i 

Let us first prove inclusion (1): 

C(cr n+ i, L n , cr n , . . . , Li, ai) C cl o conv o p(a n+1 .L* n .a n . . . L\.ai) (1) 

If n — 0, inclusion is immediate. Assume that n > 1 and let i G {1 . . .n}, we 
have just to show that JR_.r CT . oclo conv o£(L*) C cloconvop(L). Naturally, 
if = 0, this inclusion is immediate. Otherwise, let w G L*\{e}. For any 
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k G N, we have a n+ \ . . . Oi + \.w k .o~i . . . G\ C L. From the following equality, we 
get R_.r o . i ... CT1 o cl o conv o C cl o conv o p(L): 

p{a n+l . ..a w .w k .tn . . .ffi) = r CTn+ ,.. CT1 ((l - r^^.r^.^M) 

In particular, we have proved inclusion (1). Let us prove the converse inclu- 
sion. Consider a sequence (u , j)i<j< n +i such that Wi G L*\{e}. An immediate 
induction over n > 0, proves the following equality: 

/n+l 

pCffn+a.Wn+l-ffn+i . • Wl) = 1^,..^ £ r l«-+x-«*+xl . (1 - r^T"*^ O £{ Wi 

\i=l 

\ s r \w n+ i...w i+1 \ ^ _ r \wi\j e ]g>^ we (j e( j uce the following inclusion: 
cl o conv o p(a n+1 .(L* n \{e}).a n . . . (L\\{e}).ai) C C(a n+1 , L n , a nj . . . , L 1} a x ) 

Naturally, from the previous inclusions taken over n > 0, we deduce the con- 
verse inclusion of (1). □ 

The previous proposition explains why we are interested in computing cl o 
conv o where L is a regular language. In fact, we have the following 

lemma. 

Lemma 3 For any L C £* m; we have conv o = conv o £(L). 

PROOF. From L C L*, we deduce the inclusion conv o £(L) C conv o 
Let us prove the converse inclusion. Let io G L*\{e}. There exists a sequence 
Oi, 0fc of fc > 1 words in £\{e} such that w = o\ . . . 0&. An immediate 
induction over fc > 1 proves the following equality: 

k r Wi\ _ i 

e(0i...0 fc ) = E^ 1 - CT '- ll r i^i_ 1 -^) 

As Etir^-^- ^i^tu! = 1 and r ki-^-il -^g=^ G R_, we deduce that 
G conv o £(L). We deduce £(L*) C conv o and by minimality of the 
convex hull of £(£*), we get conv o C conv o □ 

Once again, we use the fact that a regular language L can be decomposed into 
a finite union of languages of the form cr n+1 .L* .a n . • • • L\.<j\. 

Proposition 4 Let us consider a language L = o n +\.L* n .o n . . . L\.o\ where 
n > 0, Oi G S* m and C S* m . We /iawe £/ie following equality: 

cl o conv o f (L) = cl o conv ^({<r n+ i . . . 0i}) (J T" 1 CTi o cl o conv o £(L;)^J 
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PROOF. Let us consider a language L of the form L = o 2 .L\.o\ where 
<Ji,a 2 G S* m and L\ C S* m and let us prove the proposition for L. Remark 
that if Li\{e\ = or if a 2 -&i = e, lemma 3 proves the proposition. So, we can 
assume that Li\{e} ^ and o 2 .0\ ^ e. Let C = cl o conv(^({o- 2 .ai}) U L^ 1 o 
cl o conv o £(Li)) and consider a G £i\{e} and fceN. We have the following 
equality: 

i(a 2 .a k = r |^| +fc , w _ -^ (^i) + rk2 , CT1 | +fc . M _ 1 -^(C^)) 

In particular we deduce that L^ 1 o £(cr) e cl o £((7 2 .(7*.(7i). From o 2 .o* .0\ C L, 
we get r^o^Li) C clo£(L). And by minimality of the closure and the convex 
hull, we get T^ 1 o cl o conv o £(Li) C cl o conv o £(L). From e 2 .oi G L, we 
deduce that £({cr2-o"i}) C We obtain C C cloconvo£(L). Let us prove the 

converse inclusion. Consider a £ L. There exists weL] such that a = a 2 .w.o\. 
If w = e then £(cr) = £(<7i.<7 2 ) G C. Otherwise, lemma 3 proves that £(u>) G 
conv o As = -gg^.fa.aj + ^iSj -^^H), we 

deduce that {(a) G C. We deduce that the other inclusion cloconvo£(L) C C. 
Therefore, the proposition is proved for L. 

Now, assume the proposition proved for an integer n > 1 and let us consider 
a language L = a n+2 .L* n+1 .a n+ i . . . L\.a x where crj G S* m and L; C £* m and 
let us prove the proposition for L. Consider w n+ i G L* l+l . As the proposition 
is proved for n, we deduce the following equality: 

cl o conv o £(a n+2 .w n+1 .a n+1 .L* n .a n . . . L\.ai) 

= cl o conv ^({a n+2 .w n+1 .a n+1 . . . <n}) JJ L^ 1 CT1 o cl o conv o f (L*)^ 

We get in particular the following equality: 
cl o conv o ^(a n+2 .L* n+1 .a n+1 .L* n .a n . . . Lj.oi) 
= cl o conv ^(a n+2 .L* n+1 .a n+1 . . . <ti}) |J r" 1 ai o cl o conv o 

= cl o conv ^cl o conv o f (a n+2 .L* +1 .<7„+i . . . a{) JJ T" 1 CTi o cl o conv o £(L;)^J 

As the proposition is proved in the case n — 1, we also get the following 
equality: 

cl o conv o £(a n+2 .L* n+1 .a n+1 . . . ai) 

= cl o conv (£({a n+2 . . . a ± }) U r^ +i CTi o cl o conv o £(L„+i)) 

The two previous equality proved the proposition for L. By induction over 
n > 1, we have proved the proposition for any n > 1. □ 
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We can now prove our main result that extends 
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Theorem 5 The convex hull of a regular set of integer vectors X = p(L) is 
polyhedral and a finite set of rays R that represents conv(X) can be computed 
in exponential time from any regular expression that defines L. 

PROOF. Let C be a polyhedral convex set represented by a finite set of 
rays R. We know that for any w G S* m , the convex sets Y W (C) and T~ 1 ((J) 
are polyhedral and respectively represented by {(r^La + c.p(w), c); (a, c) G 
R} and {(a — c.p{w), rKc); (a,c) G R}. Moreover, we also know that the 
convex set cl(M_.C) is polyhedral and represented by the finite set of rays 
{(0, 1)} U {(—a, 0); (a, c) G R}. By applying propositions 2 and 4 and lemma 
3 over a regular expression that represents a regular language L, we deduce 
that cl o conv o p[L) is polyhedral and represented by a finite set of rays 
R computable in exponential time from any regular expression that defines 
L. □ 
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